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Abstract 

We investigate the scaling properties of the two-dimensional (2D) Anderson model of localization with purely off- 
diagonal disorder (random hopping). In particular, we show that for small energies the infinite-size localization 
lengths as computed from transfer-matrix methods together with finite-size scaling diverge with a power-law 
behavior. The corresponding exponents seem to depend on the strength and the type of disorder chosen. 
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1. Introduction 

Of paramount importance for the theory of dis- 
ordered systems and the concept of Anderson lo- 
calization [1-5] is the scaling theory of localization 
as proposed in 1979 [6]. Especially in 2D, this the- 
ory predicts the absence of a disorder-driven MIT 
for generic situations such that all states remain 
localized and the system is an insulator [7-9] . How- 
ever, already early [10,11] it was suggested that 
an Anderson model of localization with purely off- 
diagonal disorder might violate this general state- 
ment since non-localized states were found at the 
band center [12-16]. Further numerical investiga- 
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tions in recent years [17-21] have uncovered ad- 
ditional evidence that the localization properties 
at E = are special. In particular, it was found 
that the divergence in the density of states DOS 
is accompanied by a divergence of the localization 
lengths A [17,18]. This divergence does not violate 
the scaling arguments [22], since it can be shown 
that its scaling properties are compatible with crit- 
ical states only [18] , i.e., there are no truly extended 
states at E = 0. Of importance for the model is a 
very special symmetry around E — which holds 
in the bipartite case of an even number of sites 
[22,23]. Then the spectrum is symmetric such that 
for every eigenenergy Ei < there is also a state 
with energy Ei > 0. This situation is connected 
with a so-called chiral universality class. Further- 
more, the model is closely connected to the random 
flux model studied in the quantum-Hall situation 
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Fig. 1. Reduced localization length A/M for various system 
sizes M of a box t distribution with c = 0. Symbols indicate 
different energies ranging from 0.025 (o), 0.0225 (□) to 
2 X 10~^ (□). The data right of the broken line have been 
used in FSS. 

where the off-diagonal disorder is due to a random 
magnetic flux through the 2D plaquettes. 

Thus although we do not have a true MIT, we 
nevertheless have a transition from localized via 
delocalized to localized behavior as we sweep the 
energy through E = 0. We consider a single elec- 
tron on the 2D lattice with N sites described by 
the Anderson Hamiltonian 



H 



N N 



(1) 



where \i) denotes the electron at site i. The onsite 
energies Cj are set to and the off-diagonal disor- 
der is introduced by choosing random hopping el- 
ements tij between nearest neighbor sites. 

We test three different distributions of Uj : 

(i) a rectangular distribution [17] 



PiUj) 



l/W if \tij - c| < w/2, 
otherwise, 
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Fig. 2. Reduced localization length A/M for various odd 
system sizes M of a box t distribution with c = 0. Symbols 
indicate different energies ranging from 0.025 (o), 0.0225 
(□) to 2 X 10~^ (□). The data right of the broken line have 
been used in FSS. 



(ii) a Gaussian distribution 
1 



PiUj) = 



: exp 



(tij - c)2 



2a2 



(iii) a rectangular distribution of the logarithm of 
tij [16] 



P{lnUj/to) 



1/w ff \lntij/to\ < w/2, 
otherwise. 



The logarithmic distribution appears more suited 
to model actual physical systems [16] . We also note 
that the logarithmic distribution avoids problems 
with zero t elements and thus there is no need to 
introduce an artificial lower cutoff as for the box 
and Gaussian distributions [17]. Furthermore, the 
box and Gaussian distributions will usually have 
negative t values which correspond to a rather ar- 
tificial phase shift. 

In the case of rectangular and normal distribu- 
tions we set the width w and the standard devia- 
tion 0- to 1 and change the center c of the distribu- 
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Fig. 3. Reduced localization length \/M for various system 
sizes M of a Gaussian f distribution with c = 0.25. Symbols 
indicate different energies ranging from 0.03 (•), 0.0275 
(■) to 2 X f0~^ (a). The data right of the broken line have 
been used in FSS. 

tion. In the case of the logarithmic t distribution 
to = 1 sets the energy scale and we change the 
disorder width w. Values of the parameters were 
c = 0, 0.25, 0.5, and 1.0 for the rectangular distri- 
bution; c = and c = 0.25 for the Gaussian distri- 
bution and w = 2, 6, and 10 for the logarithmic t 
distribution. 



2. Computation of the localization lengths 

The transfer-matrix method [24,25] was used to 
compute the localization lengths for strips of vari- 
ous widths M up to M = 100 in the energy interval 
2 X 10~5 <E< 0.2048. In Figs. 1, 2, 3, and 4 we 
show the system size dependence for, e.g., special 
values of c and w and all three disorder distribu- 
tions. The accuracy of our results is 0.1 — 0.3% or 
1% depending on the disorder distribution and the 
values of parameters, see Table 1 for actual param- 
eter values. 




Fig. 4. Reduced localization length A/M for various system 
sizes M of a logarithmic t distribution with w = 2. Symbols 
indicate different energies ranging from 0.2048 (o), 0.1024 
(□) to 2 X 10-" (o). 

Next, the finite-size-scaling analysis of Ref. 
[25] was applied to the data. The calculated lo- 
calization lengths usually increase as the energy 
approaches 0. Only, for small even width values 
(10,20) it decreases significantly close to E = Q 
[21] which makes finite-size scaling impossible. 
Therefore the smallest system sizes were dropped 
during the finite-size scaling procedure. Results 
for the finite-size scaling curves are shown in Fig. 
5 for the three different distributions. 



3. Critical exponents 

One expects that the scaling parameters ^ ob- 
tained from finite-size scaling diverge close to E = 
[26]. However, the precise functional form of this 
divergence is not yet know. In Ref. [26] it has been 
suggested that for energies E > E* the divergence 
can be described by a power law as 
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Fig. 5. Finitc-sizc-scaling plots for box (c = 0, M € [30, 80], 
large open symbols), Gaussian (c = 0, M G [20,60], filled 
symbols), and logarithmic {w = 2, M E [30,100], small 
open symbols) t distributions. 

with tlic critical exponent v. For even smaller 
\E\ <C E* , this behavior should then change to 



^{E) oc exp 



InEo/E 



A 



(3) 



with constants Eq and A given by the renormal- 
ization group flow [26]. Double-logarithmic plots 

of ^ vs. E in Figs. 6, 7 and 8 confirm the power-law 
behavior with reasonable accuracy down to i? « 
10~*. For smaller values it has been shown already 
in Ref. [21] that a new behavior is to be expected. 

Table 1 collects the vahies of the critical expo- 
nent obtained for different disorders. In the case 
of the logarithmic t distribution and w = 10 the 
power-law divergence fails, therefore the expo- 
nent was not calculated. From Table 1, it can be 
easily seen that all calculated values are in the 
range 0.2 < < 0.5. The exponent is apparently 
not universal but seems to depend on the kind 
of disorder and the actual value of parameters; 
for stronger disorders v becomes smaller (for the 
logarithmic t distribution the disorder strength 
increases with w [16], for the rectangular distri- 
bution the strongest disorder appears at c = 0.25 
[17]). This non- universality is in agreement with 
the results of Ref. [26]. 



3 - 



O 




-5 



o 


C=0 


+ 


c=0, odd M 


□ 


c=0.25 




c=0.5 


A 


c=1 


1 1 


-4 



Fig. 6. Variation of the infinite-size localization length f 
with E for box distributions. The inset shows the t distri- 
bution for c = 0. 




Fig. 7. Variation of the infinite-size localization length { 
with E for Gaussian distributions. The inset shows the t 
distribution for c = 0. 

As the localization lengths calculated for odd 
and even strip widths may exhibit different be- 
havior [16,21] we repeated the calculations also for 
odd-width systems for chosen parameter sets for 
all distributions. Fig. 2 shows an example of re- 
duced localization lenghts for rectangular distribu- 
tion. In contrast to the even-width systems (Fig. 
1) the localization lengths do not decrease signif- 
icantly close to the = 0. We have attributed 
the difii'erent behaviors for odd and even system 
sizes to different structures in the density of states 
[21,27]. Nevertheless, the exponent in both cases is 
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Table 1 



Estimated values of the exponents 


of the localization lengths 


for various disorder strengths and distributions. The error 


bars represent the standard deviations from the power-law fit and should be increased by at least 
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for a reliable representation of the actual errors. 






disorder 


parameters 


accuracy 


sizes used 
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distribution 
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in finite-size scaling 


exponent 
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Fig. 8. Variation of the infinite-size localization length ^ 
with E for logarithmic distributions. The inset shows the 
t distribution for ui = 2. 

within the error bars the same (cp. Table 1). This 
is also true for Gaussian and logarithmic disorder 

distributions (sec Tabic 1). therefore, at least for 
the investigated disorder strengths, the difference 
in exponents of localization lengths is negligible. 



diagonal Anderson model of localization in 2D can 
be described by a set of exponents that model the 
divergence of the localization lengths ^ at = 0. 
Note that these exponents are in reasonable agree- 
ment with the exponent 0.5 first estimated for the 
scaling of the participation numbers in Ref. [17]. 
Down to « 2 X 10~^ in Figs. 6, 7 the power- 
law behavior can model the data reasonably well. 
Thus we expect the crossover predicted in Ref. 
[26] to appear at smaller energies. We find that the 
exponents depend on the strength and distribu- 
tion of the off-diagonal disorder also in agreement 
with Ref. [26]. Currently, we are extending these 
calculations to smaller energies. 

We note that it might be interesting to also in- 
vestigate the situation in honeycomb lattices [28] , 
where the van Hove singularity of the square lat- 
tice &tE = does not interfere with the divergence 
due to the bipartiteness which is of interest here. 
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4. Conclusions 

Our results suggest that the localization- 
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